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Abstract 

We study the stability and the modes of non - isothermal coronal loop models with different 
intensity values of the equilibrium twisted magnetic field. We use an energy principle obtained via 
non - equilibrium thermodynamic arguments. The principle is expressed in terms of Hermitian 
operators and allows to consider together the coupled system of equations: the balance of energy 
equation and the equation of motion, to obtain modes and eigenmodes in a spectrum ranging from 
short to long-wavelength disturbances without having to use weak varying approximations of the 



C"^^ '_ equilibrium parameters. Long-wavelength perturbations introduce additional difficulties because 



the inhomogeneous nature of the medium determines disturbances leading to continuous intervals 

of eigenfrequencies which cannot be considered as purely sinusoidal. We analyze the modification of 

• 1—1 . 

r> I periods, modes structure and stability when the helicity, the magnetic field strength and the radius 

S ' 

of the fluxtube are varied. The efficiency of the damping due to the resonant absorption mechanism 

is analyzed in a context of modes that can either impulsively release or storage magnetic energy. We 

find that the onset of the instability is associated to a critical value of the helicity and that the 

magnetic energy content has a determinant role on the instability of the system with respect to 

the stabilizing effect of the resonant absorption mechanism. 
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I. INTRODUCTION 

A. Variational principle 

A crucial requirement for any theoretical model of coronal structures is to give account 
of the stability and evolution of far-from-equilibrium states which are responsible of the 
characteristic rich topology and dynamics of the solar corona. This implies to consider the 
coupling of thermal and mechanical equations. Different stability analysis of solar structures 
can be found in the literature, generally restricted to special types of perturbations and 
specific equilibrium models. These includes, models that consider adiabatic configuration 



Q 



such as the ones analyzed via the classical criterion of Bernstein j^ or those that presuppose 
static equilibrium and analyze thermal stability. In the application of Bernstein's criterion, 
the adiabatic assumption implies that the energy balance equation is not required and 
thus dissipation is impossible. Also the assumption of static models is a strong, and often 
unjustified, restriction for open systems. 

In this paper we apply an energy principle to analyze the stability of solar coronal loops 
when helical modes are present. The principle was obtained in previous papers (Paper I: |5|; 



6j; see also 20|) using a general procedure of irreversible thermodjTiamics -based on firmly 
established thermodynamic laws- that can be understood as an extension of Bernstein's 
MHD principle to situations far from thermodynamic equilibrium. 

In Paper I and in |6!] we showed how to obtain the variational principle for solar coronal 
structures from the equations that describe the dynamics of the system. The method consists 
of obtaining a Lyapunov function, also known as generalized potential, that represents the 
mathematical expression of the stability conditions. The principle is subject to physically 
reasonable requirements of hermiticity and antihermiticity over the matrices. For a more 
detailed presentation see Paper I and the references therein. 

B. Solar coronal loops 

MHD loop oscillations in the corona are known to be strongly damped, mostly having 
decaying times of few periods Np ^ 2 — 7 periods. While thermal conduction, with the 
contribution of radiative cooling mechanisms, could be the main cause of the damping of 
pure MHD slow magnetoacoustic mode oscillations they are unimportant for the MHD fast 



modes. Resonant absorption and phase mixing seem more promising in giving account of 



the rapid decay (111]; hereafter HG, [10]) of the ideal fast oscillations of these strongly 
inhomogeneous and structured plasma systems. Inhomogeneous equilibrium distributions 
of plasma density and temperature varying continuously across the magnetic field led to 
plasma waves with continuous intervals of eigenfrequencies. The occurrence of the Alfven 
ideal MHD continuum in a thin edge layer is derived from the highly anisotropic charac- 
ter of the fast magnetoacoustic waves giving rise to a peak of the amplitudes where the 
perturbation develops large gradients and the absorption has maxims. However, there is 
another type of continuum commonly known as slow magnetosonic continuum associated to 
the inhomogeneity of the equilibrium parameters along the axis of the loop (see Paper I). 
This inhomogeneities are associated, for example, to changes in the density concentration 
at the loop basis. If the magnetic field is twisted the inhomogeneities led to the coupling of 
Alfven and slow magnetosonic continuum modes ((3|). 

The resonant absorption mechanism of wave heating consists on the non-dissipative trans- 
ference of wave energy from the collective line-tied wave with fast discrete eigenvalues (kinetic 
energy of the fast radial component) to a local resonant mode in the Alfven continuum, (ki- 
netic energy of the azimuthal component), which is then dissipated in an enhanced manner. 
Then, the continuum oscillations are converted into heat by dissipative processes; as the 
medium has large gradients in the Alfven speed, the oscillations of neighboring field lines 
become out of phase and shear Alfven waves lead to enhanced viscous and ohmic dissipa- 
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tion (see [12| for the linear regime and [iJ] for the nonlinear one). The mode conversion 
from the collective to the local mode occurs in a time that is non-dissipative and generally 
much shorter than the second time scale which is related to the dissipative damping of the 



small-scale perturbations of the local mode in the resonance layer ( 18|]; 22|). 

The whole temporal pattern description of modes that exhibit a combination of global 
(discrete line-tied fast eigenmode) and localize (Alfven continuum mode) behaviour is known 
as quasi-mode. Moreover, the mixed nature of the modes is not only due to the temporal 
behaviour but also to the boundary value problem giving rise to a spatial behaviour which 
is also of a mixed nature, i.e. coronal loops with line-tying constraints cannot support 
pure waves: Alfven, slow or fast magnetoacoustic modes. HG studied the mixed spectral 
description of coronal loops (i.e. the resulting superposition of basic waves which adjust 
the line-tied condition) without assuming a straight magnetic field and forcing the loop 



to follow the photosplieric velocity perturbations. They found that pure Alfven and pure 
slow modes are obtained as singular limiting cases of cluster spectra of Alfven-fast or slow- 
fast modes, where the fast components are localized in a photospheric boundary associated 
to the line-tied condition: the coronal part of the loop acting as a resonant cavity of large 
Alfven components and fast components, with a small but rapidly varying amplitude, located 
in the photospheric boundary layer. They found that heating of coronal loops by resonant 
absorption is due to the line-tied Alfven continuum which no longer depends on the poloidal 
magnetic field and that the corresponding eigenmodes have a global ballooning feature which 
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is characterized by an accumulation point given by the Alfven frequency. In [9| (hereafter 
GH), a variational principle, based in Bernstein's principle, was obtained to derive the 
Alfven and slow continuum frequencies in a line-tied inhomogeneous cylinder. Stability 
considerations led them to conclude the global stability of coronal loops. 

In this paper, following results of Paper I we apply our energy principle to consider the 
stability and mode structure of loop inhomogeneous coronal models with non-vanishing 
helicity. Our principle has the advantages that it does not require a WKB approximation 
and that, as was mentioned, it allows the consideration of the coupling of the thermal and 
mechanical equations that are necessary to analyze far from equilibrium states. 

II. THE MHD STABILITY CRITERION FOR CORONAL STRUCTURES 

Solar coronal conditions with large Reynolds numbers are well fitted by ideal MHD plasma 
models (i.e. infinite electrical conductivity a ^ 1 leading to vanishing viscosity and ohmic 
dissipation). Thus, the fundamental equations considered are the mass conservation equa- 
tion, the perfect gas law or state equation for a fully ionized H plasma and the induction 
equation, with vanishing magnetic diffusivity due to the conductivity properties. The energy 
balance equation takes the form: 

P' ^ —(^) = -V -F.-U + H (1) 

Fc is the heat flux due to particle conduction along the loop, Lr is the net radiation flux 
and H the heating function which was chosen as in Paper I: H = hp + Hq. Eq. [T] expresses 
the fact that the gain in particle energy (internal plus kinetic) is due to the external heating 
sources represented by the heating function, heat flow and radiation losses; all other heating 



sources were considered as vanishing terms implying that the optically thin assumption 
holds. Note that the non-ideal contribution in the energy equation (L) is associated to the 
open character of the loop system. 

Once the linearization around a nonlinear equilibrium or stationary state is performed, 
and after a straightforward manipulation procedure where the hermiticity requirements are 
fulfilled the generalized energy principle and the respective frequencies are obtained (Paper 
I and [6]) as: 

5'^Wp = - /(f /9Ff + T^ATi + T*Bi- CBTi)d^x > 0. (2) 

2 _ Ii^*PF^+ T*ATi + T*Bi- i*BT*)d^x 
u) — — — (3) 

with the same normalization condition as in Paper I. F is the known Bernstein operator 
for the system, ^ and Ti are the motion and temperature perturbations and operators A 
and B are as in Paper I. For the non-dissipative cases (L = or equivalently Ti = 0), last 
expressions (discarding the presence of factor (3 which appears in the equations to fit the 
Hermitian and anti-Hermitian conditions) are reduced to the well-known Bernstein MHD 
energy principle and its respective frequencies. 

III. APPLICATION TO AN INHOMOGENEOUS LOOP MODEL WITH NON- 
VANISHING HELICITY 

On one hand, the azimuthal component of the loop perturbation is believed to be one 
of the principle responsible of resonant absorption and damping of ideal oscillations; on 
the other, this component is associated to the storage of magnetic energy in systems with 
non-vanishing helicity which eventually is released by instabilities. Thus, we are interested 
in analyzing the changes produced in the stability of non-homogeneous loops subject to 
helical perturbations. This is, loops with inhomogeneous distributions of plasma density 
and temperatures subject to body modes and with non-vanishing helicity. In this case, the 
Alfven, slow and fast magnetoacoustic cylinder modes cannot longer be associated to the 
azimuthal, longitudinal and radial components respectively. The observational importance 
of helical modes cannot be neglected and it is poorly known how helicity affects important 
physical features of mode oscillations (e.g., damping mechanisms, stability and periods). 
However, a mode classification can be accomplished via the analysis of the mode variations. 



described in an orthogonal basis, while helicity is varied. The basis is formed by the orthog- 
onal displacements: parallel and perpendicular to the magnetic field and the radial (and 
perpendicular to the surface of the tube) one, of observational interest. 

The fundamental modes are generally observationally and energetically more important 
than their harmonics. For these global modes the inhomogeneous nature of the medium 
cannot be ignored and it determines the structure of the disturbance which cannot be taken 
as sinusoidal, making the traditional normal mode analysis useless for this treatment (sinu- 
soidal dependence with constant coefficients), i.e. at least a WKB approximation, of weakly 
varying parameters compared to a typical wavelength, is required. Moreover, the occurrence 
of either an infinitely degenerate eigenvalue or an accumulation point giving rise to a contin- 
uous spectrum are associated to inhomogeneities. We consider two types of inhomogeneities: 
the inhomogeneity of the equilibrium parameters along the loop axis, and the inhomogeneity 
across the loop axis when the radius is varied. As a first order approximation we neglect 
the effect of gravitational stratification and thus confine the analysis to characteristic spa- 
tial scales lower than the pressure scale height in the solar corona. In order to analyze the 
stabihty and to obtain the frequencies and modes the physical quantities in eq. [2] and eq. H] 
must be calculated along the loop structure. 

A. Mechanical equilibrium 

To determine an equilibrium configuration we assume force-free equations. This as- 
sumption is justified for coronal conditions due to the fact that in plasmas with low (3 (gas 
pressure over the magnetic pressure) the pressure gradient can be neglected in comparison to 



the Lorentz force. For the chromosphere and the photosphere the 



brce-free approximation 



may not be a good one. However, it is a widespread supposition [l9|: perturbed systems are 
believed to relax to new force-free, minimum energy states and chromospheric conditions 
seem to be well fitted to force-free models from 4. 10^ m [2| (Chapter 5). 

Coronal loops are generally modeled as thin cylindrical fiuxtubes where the curvature and 
related forces can be neglected so the cylindrical geometry can be applied. The fiuxtube is 
assumed as line-tied to the photospheric plasma through its footpoints which are forced to 
follow the photospheric velocity perturbations. The random velocity field creates vorticity 
generally twisting the coronal fiuxtubes. Thus, a relation between the helical twist and 



the force-free parameter can be derived as follows (e.g. [2l|]). The coronal loop model is 
obtained from the equations 



VxBo = a(r)Bo j x Bo = 0. 



(4) 



Also, since Bq is force-free, Vp = everywhere and thus has a constant value along the loop. 
We consider a straight cylinder with a nonuniform distribution of density and temperature 
and a resulting uniform twist over an initially non-rotated field B = (0, 0, B^) yielding the 
unperturbed magnetic field 



Bo = (Br, B^, Bz) = Bq{0. 



br 1 , 



with A = 1 + 6^r^ and b = 2'U.Nt/ L {Nt number of turns over the cylinder length L). Then, 

(5) 



Biji 



dz 



r27iNt , , , 26 
= brair) = — - 



We assume a given value of the cylinder radius r = R, thus the line element results a function 
of the coordinate z: s = s{z). The dependence with the radial component will be taken into 
account by considering different values of the radius R. 



ds^ = R^d(f)'^ + dz^ = (l + R^b"^) dz^ = Adz^ 



(6) 



B. Thermal equilibrium 



The thermal equilibrium is obtained, as in Paper I, assuming L = in the balance 
energy equation (eq. [1]) . The procedure developed consists in obtaining the function of 
the temperature along the arc element s by integrating eq. [H with the constraint L = and 
replacing border conditions: the temperature at the bottom T^ = lO'^K and the temperature 
at the top Tj = 10^ K. The known expression (see chapter 6 of Priest 15|) is obtained 



dT 
ds 



P^X 



2klko(a + I 



-T"-! 



which has to be inverted to obtain T = f~^{s) [1] as 






(7) 



^'^ - A 

ds 



dMy dv 
dv dT 



where B^(^,g)= f t^-^ {I - tf-^ dt 
2 Jo 



(8) 



with 



p = ^; ^ = 1 - (^r-; g = (f + ^)(2 - «) + 1 



A={2- a)Tr^ {{p\)/{2ko{a + ^)4))^. 



We use a = — |sog = |to numerically calculate the modes, 

. = ^B.(1,|)^^=|T3(^)V2. 

Also, from boundary conditions v = 0, thus the constant value of the heating function results 

Ho = 7p\TrV {8kl{a + §)) . 

C. The perturbation 

To calculate the stability and the structure of the modes the general perturbation along 
the equilibrium magnetic field is written 

f = [C.(r, z)et + za(r, z)e^ + C.(r, z)e,y"''^ T, = T^{r, z)e'^'^ (9) 

with r = R. The dependence only appears in the exponents that multiply the perturbation; 
the integration with respect to this coordinate is straightforward. Then, representing the 
equilibrium functions of the different quantities with a sub-index, defining 
et = {Rbe^ + e,)/VA Vy = et{et ■ V) pt = ^ 

with e^, Bz the cylindrical versors and e^ the tangential versor, we obtain a non-dimensional 
expression for the energy principle of eq. EJ 

S^Wp = S^Wc + 6^Wm + S^Whc + S^Wr (10) 

where S'^Wc is the generalized potential energy associated to compressional terms, S^Wm 
corresponds to the magnetic contributions, 5'^Whc corresponds to the heat conduction terms 
and 5'^Wr to the radiative contributions. The explicit form of these functions are given in 
the Appendix. The Bernstein's generalized potential energy corresponds to the magnetic 
contribution and part of the compressional one. In the generalized version of the energy 
principle additional terms appear in the 6'^Wc term and also S'^Whc and 6'^Wr are entirely 
new terms. 
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IV. RESULTS AND DISCUSSION 

Convective motion of the photosphere is beheved to provide the energy that is storage 
in twisted magnetic coronal fields allowing the presence of long-lived coronal structures 



until it is released by instabilities ( 16j; 23|). On the other hand, continuous spectra are 



generally associated to stability. An accepted conjecture establishes that unstable modes 
have a discrete spectrum (see 7| or [l5|). There are two types of possible continuous spectra 
in this problem. The inhomogeneous character of the equilibrium parameters along the 
loop axis can lead to a continuum that couples to the Alfven continuum ^] e.g, when 
the disturbances considered are comparable to the inhomogeneous characteristic wavelength 
stable eigenvalues can give rise to a continuous spectrum {L/2, the equilibrium structure in 
the z component). This is the case studied in Paper I. On the other hand, GH established, 
for non vanishing helicity systems, that there is a continuous spectrum associated with the 
line-tied Alfven resonance leading to the damping and heating by the resonant absorption 
mechanism and thus, directly relate to the stability of loops. They also pointed out how to 
obtain the resonant singular limit uj, from the class of physically permissible solutions, 

uj{r) = -^ — . (11) 

This resonance results because of the absence of an explicit dependence on the azimuthal 
magnetic field component {B^). 

Thus, in order to understand in which conditions which mechanism can dominate and 
give account of the different scenarios i.e., the driving of the instability or the damping 
of mode oscillations, it is critical to gain knowledge about the dynamics and energetic 
contribution of twisted structures. Yet, the implications of the twisting in theoretical and 
observational descriptions are poorly known; e.g., there is no clarity about the modification 
of the dispersion relation and observational data are indirectly inferred. 

In this paper we focused our attention to describe the changes in periods, stability and 
mode structure of coronal loops when the helicity, the magnetic field intensity and the radius 
are varied. For loops with vanishing helicity it is well established that the Alfven line-tied 
resonance continuum is responsible of the damping of kink {171 = 1) quasi-modes via the 
transfer of energy from the radial component into the azimuthal one, i.e., from discrete 
global modes into the local continuum modes where phase mixing can take place. Still, 



the twisting of the magnetic field leads to the coupling of MHD cylindrical modes making 
difficult to provide a classification in terms of the behavior of pure-like modes. 

In order to calculate modes and frequencies we followed the schematic procedure de- 

n 

scribed in Paper 1 and in [8|. We used a symbolic manipulation program to integrate the 
equations. S'^Wp and the perturbations were expanded in a six dimensional-Fourier basis 
on the independent coordinate z that adjusts to border conditions, i.e., the four pertur- 
bated components (eq. [9]) were expanded in a six mode basis to obtain 24 eigenvalues and 
eigenvectors for each of the helicity and the magnetic field values. Only the first eighteen 
eigenvalues were considered (the others are more than two order of magnitude smaller and 
accumulate at zero; the eigenvectors are also vanishingly small). Thus, a quadratic form 
for 6'^Wp was obtained and was minimized with the Ritz variational procedure. A matrix 
discrete eigenvalue problem subject to a normalization constraint was obtained. From the 
resulting modes and the generalized potential energy (eq. [2]): 6^Wp > the stabihty of each 
mode was determined. 

The coronal loop parameters used were: L = 10^^ cm (or L = lOOMm), Tb = W^K 
Tt = lO^K He = 10^cm~^ electron number density pt = 2kBTt ; pt = rUpPt/kBTt. Frequen- 
cies and modes were calculated for two different values of the magnetic field: Bq = lOG and 
Bo = lOOG, and for three different values of the helicity 6 = (3.1 10"^; 3.1 10~^; 1.9 10"^) 
which correspond to the adimentional values: ba = (2.8; 28; 170) with A^^ = (0.45; 4.3; 13.7), 
Nt : the number of turns over the cylinder length. These helicity values defined as weak, mod- 
erated and strong helicity respectively correspond to the classification given in [2] (Chapter 
5). The adimentional radius was initially chosen a.s R = 0.01. In what follows we summarize 
the conclusions obtained from the data analysis which are displayed in three tables. 

Table 1 shows the periods (in minutes) for weak, moderate and strong helicity for two 
values of the magnetic field intensity {Bq = lOG and Bq = lOOG (left and right panel 
respectively). S and U letters indicate the stable-unstable character of the modes. From 
the table we see that: 

I) Weak helicity modes are stable. This is in accordance with the analytic results by Rud- 
erman [l9| who studied nonaxisymmetric oscillations of a thin twisted magnetic tube with 
fixed ends in a zero-beta plasma. 

II) Higher modes have an accumulation point at zero, indicating the presence of a continuum 
spectra of stable modes (as in Paper I). Note that, calculus performed via discrete basis, 
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as in our case, give spectra that are necessarily discrete. Thus, an accumulation of discrete 
eigenvalues suggests a stable continuum spectrum. 

III) The Bq = lOG case has larger periods than the Bq = lOOG one. For moderate and strong 
helicity the eigenvalues of the first panel follow a scaling law with that of the second one i.e., 
they scale with the magnetic field intensity exactly as the Alfven speed does Piqg — 10-PiooG- 

IV) As HG and GH, we note a clustering of the spectra associated to the change from real to 
imaginary eigenvalues (and viceversa). There is a pronounced change (in the spacing of the 
periods or/and in the stability) from the sixth mode to the seventh mode. This is noted by a 
double line in Table 1 and related to the importance of the parallel component with respect 
to the perpendicular component (see Table 2). Up to period number ten real - imaginary 
eigenvalues of the first panel {Bq = lOG) correspond to real-imaginary ones of the second 
panel [Bq = lOOG). Also, excepting large order periods n > 10, when the helicity is increased 
from weak to moderate the imaginary stable eigenvalues turn to imaginary unstable ones. 
For Bq = lOG and weak and moderate helicity cases there are five different groups of periods 
(Pi — Pq] P-j — Pio; Pii — P12; Pi3 — P14; Pi5 — Pis) (see also Table 2). The clustering is more 
difficult to establish i.e., the differences are less pronounced, with increasing magnetic field 
intensity and larger order periods. 

In order to compare our results with those given by these authors we calculated the expres- 
sion eq. [11] for our modes. We found that, all periods excepting Pi — Pq weak helicity modes 
satisfy the relation and thus, they belong to the Alfven continuum spectrum justifying the 
scaling law described in III. As HG, we conclude that the change in the real-complex char- 
acter of the Pe — P-j eigenvalues is associated to the existence of an accumulation point of the 
resonant Alfven continuum, however we find that this change is not necessarily related to a 
change in the stability as they claimed. Note that all modes with weak helicity are stable 
(even the imaginary ones); in all the other cases the imaginary character of the eigenvalues 
is associated to instability. Yet, the continuum stable eigenvalue conjecture is here still valid 







151]; it applies to a spectrum with an accumulation point in zero; we found stable modes 
for all the helicity values and for the two magnetic field values with Pn>i4. Note that the 
analysis of stable modes is still of interest because depending on the relative characteris- 
tic times of stable and unstable modes the stable ones could be active and accessible to 
observations. 
The presence of at least one unstable mode means that the equilibrium state is unstable. 
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Thus, taking into account the whole range of stable modes, we confirm previous results 
leading to conclude that field configurations with some degree of twisting give a stabilizing 
effect allowing the storage of magnetic energy 17|, i.e., when the helicity is augmented the 



stable weak case turns to an unstable one suggesting a critical value. 

In Paper I, we obtained only one unstable mode classified as slow magnetoacoustic mode 
due to the almost longitudinal character (parallel to the magnetic field) of the wavevector 
perturbation and to the fact that the period did not changed with the intensity of the 
magnetic field, resembling acoustic waves with sound speeds, Vg, independent of the magnetic 
field. The characteristic unstable time obtained in Paper I was r^ = 36 min, corresponding 
to a typical slow magnetoacoustic fundamental period with a characteristic wavelength of the 
order of the loop length L/2. Also, we obtained a continuous set of stable modes classified as 
fast magnetoacoustic modes due to their large value component orthogonal to the magnetic 
field and to the fact that the eigenvalues scale with the intensity of the magnetic field as 
PiiQ ~ IOPioog; thus resembling the dependence of the Alfven waves va ^ Bq. 

Table 2 (First Panel) displays the resulting features associated to the relative intensity 
of the parallel and perpendicular to the field components ((^y, ^_l, Cr) is an orthogonal basis) 
and their classification as slow-like (S) or fast-like (F). The relative phase between the 
components is also indicated in the table by P (in phase) and IP (inverted phase). Table 2 
(Second Panel) also shows the intensity relationship between the cylindrical components. In 
order to classify the modes and to compare with the slow and fast magnetoacoustic modes 
obtained in Paper I, we calculated the cylindrical mode components and also the tangential 
and normal to the field components (^y = (-Rfe^,/, + ^2)/A; ^_l = (,^(^ — RbC^z)/A). Our interest 
in the ^y, ^z, ^r, ^± and ^^ components resides in that: First, when the helicity is weak, 
the ^11 component is expected to play the slow-mode role of C,z in Paper I. Second, the ^r 
component is related to the fast modes and determines the resonant absorption mechanism 
when uniform cylindrical fiux tubes are considered by the transferring of energy to the ^^ 
component. When helicity and inhomogeneous distribution of equilibrium parameters are 
present it is worth investigating the transferring of energy from the C,r component to the 
others. In this case the resonant damping of global oscillations will occur by conversion of 
kinetic energy of the radial component into kinetic energy of the C,\\ and C,± components; 
both components forming the plane orthogonal to C,r, and equal to the plane formed by ^^ 
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and ^z- 

From the analysis of the amphtude of the components of the Pi — Pg modes with respect 
to the P7 — Pis ones in the weak hehcity case i.e., real and imaginary eigenvector respectively, 
we could classify the first ones as slow-like modes because: I) their tangential components S,\\ 
are at least an order of magnitude larger than the normal ones ^±\ II) as the helicity is weak 
^11 K, ^z and C,r ~^ O;^0 ^ 0, the wavevector is almost tangential to the magnetic field; III) 
they have a larger characteristic time and a shorter characteristic speed than the imaginary 
eigenvectors. On the contrary, imaginary eigenvalues are associated to large values of the ^r 
component and ^_l component (due to large values of ^^ (see Table 2 Second Panel)) , and 
small values of the C,\\ and C,z components. As in Paper I, when the eigenvalues change form 
real to imaginary the period strongly diminishes and a change in the type of mode from 
the slow to fast magnetoacoustic type occurs. In opposition to Paper I where the acoustic 
mode has the same eigenvalue for both magnetic field intensities, here the modes are affected 
by the strengthening of the magnetic field leading to an-order-of-magnitude shorter period 
than in the non-helicity case. The C,\\ and ^_l components are in an inverted phase for real 
eigenvector modes and in phase for imaginary eigenvector modes. 

For moderate helicity the overall description is similar but all the cases having non vanish- 
ing ^0 component and all the periods in the resonant line-tied continuum. As was mentioned, 
real-imaginary eigenvalues correspond to stable-unstable behavior. 

In the strong helicity case, as the weak and moderate ones, we note for Pi — Pg larger, 
but comparable, values of the ^y component with respect to the ^_l component. In this 
case the two components of the mode are in phase. This relationship between the C,\\ and 
^± components of Table 2 (FP), and their associated phases is found again in the modes 
with Pi5 — Pig. In spite that these features are associated to the slow magnetoacoustic 
characterization. Table 2 (SP) shows that as C,z is vanishingly small, the strong helicity case 
cannot be classify as a slow mode. 

When helicity is present the mixed character of the modes manifests itself making difficult 
to identify the components that are involved in the damping mechanism. However, taking 
into account the resonant frequency of eq. [HI we noted that (HG) all the modes, except 
those with Pi — Pg periods of the weak helicity case, have resonant frequencies suggesting 
that resonant absorption in helical modes is associated to modes with significant values of 
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^_L component. If this argument is correct we can affirm that the damping mechanism of 
body hehcal modes is associated to the transfer of kinetic energy of the radial component 
into kinetic energy of the C,± component which is not only related to the ^^ cylindrical 
contribution but also to the C,z one by the expression ,^_|_ = {C,^ — Rb^z)^- 

We also analyzed the change of the period as a function of the radius for different values of 
the helicity We found that, for weak helicity, the increasing of the radius leads to a decrease 
of the periods. This is in accordance with observations, e.g., observed sausage modes are 
associated with thicker and denser loop structures and lower periods; while in other case 
(unstable cases) the increasing of the radius leads to an increase of the period. 

Table 3 -First and Second Panel- shows the variation of the radius R with the twist bR for 

n 

weak and moderate helicity respectively. [19] has conjectured that the line-tying condition 
at the tube ends should stabilize the tube and has suggested a critical value (~ Lb < q; with 
q a positive constant and L the loop length) for the onset of instability. Also, [l3| found 
that when the helicity grows beyond a critical value, the kink isolated twisted magnetic flux 
tubes below the photosphere become unstable. In fact. Table 3 can be seen as the variation 
of R with the twist value: bR, for constant values of the helicity b in the two cases: weak 
and moderate respectively. Stability is guaranteed when the loop radius is varied between 
R = 0.01 and -R = 0.1 and the helicity is weak b = 0.05 (for almost the same value of the 
length of the loop, L). However, when the helicity is incremented to 6 = 0.5 even for the 
radius oi R = 0.01 the loop structure is unstable, thus, instability can be associated with 
the presence of helicity values larger than a critical one. 

Figure [1] shows the general potential energy for Pg and P7 in the weak and moderate 
cases. Note the change of this function when the system turns from stable to unstable, as 
hehcity is augmented i.e., from 6'^Wp > to 6^Wp < 0. Figure [1^ and Figure [lb display 
the total energy composed by the compressional, radiative, thermal and magnetic energy 
contributions of Pg mode in the weak and moderate case respectively. The same features 
but for the P7 mode are shown in Figure [1^ and Figure [1^. Figure [TJj and Figure [T]i show 
the magnetic energy content alone for Pg mode in the weak and moderate case respectively. 
Figure [Uf and Figure [Di show the magnetic energy content for P7 mode and for the weak 
and moderate case respectively. It can be seen, in this and in all the other cases, that the 



14 



magnetic energy content has a determinant role on the stabihty-instabihty of the system, 
i.e., the stabihty changes when the magnetic generahzed potential energy changes sign. 
Thus, a result of this analysis is that the stability of twisted coronal loops is fundamentally 
determined by the storing of magnetic energy, being the other contributions less significant. 
Meanwhile, when the helicity is weak or vanishingly small and the magnetic contribution has 
a stabilizing effect the other non-dominant contributions, as the non-adiabatic ones, can play 
an important role. This makes possible, for example, the damping of fast excitations due 
to resonant absorption. Yet, even when one of these contributions is unstable, stable modes 
could be active for a while if their characteristic periods are shorter than the characteristic 
time of the instability. This is the case of Paper I, where we obtained a slow mode with an 
unstable characteristic times of r ~ 36 min coexisting with stable fast modes with periods 
about P ~ 1 min; moreover, we showed that the instability can be nonlinearly saturated 
giving rise to a limit-cycle solutions, i.e., an oscillation between parallel plasma kinetic energy 
and plasma internal energy where the magnetic energy plays no relevant role. Thus, the 
contribution to the stable-unstable character of the modes is mostly due to the magnetic 
energy content and not to other energetic contributions. Note that as the balance energy 
equation takes into account non-adiabatic contributions, i.e., radiation, heat flow and heat 
function (with L = at the equilibrium), the resulting perturbations are not constrained 
to the force-free condition. So, one result of the analysis is that the pertubation energy 
contribution is mainly due to magnetic forces. Thus, for these type of twisted magnetic field 
models, non-adiabatic perturbations (e.g. thermal perturbations) and resonant absorption 
seem unimportant to guarantee stability; a loop system with weak storage of magnetic energy 
(low values of the helicity) could be released if the helicity is suddenly increased, e.g., by 
footpoint motions. Meanwhile all the "zoo" of the coronal seismology can be active and 
accessible to observations. 

V. APPENDIX: GENERALIZED POTENTIAL ENERGY TERMS 

From the procedure described above and extensively exemplified in Paper I we can obtain 
-laboriously but in a straightforward way- the explicit terms for the energy principle given 
in eq. [101 
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FIG. 1: Energy content of the sixth and seventh mode for Bq = lOG. a) Total potential energy 
and b) magnetic potential energy respectively for the sixth mode Pq = 1.23 min and for weak 
helicity. c) Total potential energy and d) magnetic potential energy respectively for the sixth 
mode Pq = 1.23 min and for moderate helicity. e) Total potential energy and f) magnetic potential 
energy respectively for the seventh mode Pj = 0.07 min and for weak helicity. g) Total potential 
energy and h) magnetic potential energy respectively for the sixth mode Pj = 0.07 min and for 
moderate helicity. 

where the right side of the equation corresponds to the compressional, magnetic, heat 
conduction and radiative contributions respectively. The compressional term 6'^Wc = 
6'^Wci + S'^Wc2 has an additional contribution {6^Wc2) with respect to Bernsteins princi- 
ple: 
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TABLE I: Eighteen first periods associated to stable (S) and unstable (U) eigenvalues (minutes) 
for A) Left panel: Bq = lOG with Al) left column: weak helicity, A2) middle column: moderate 
helicity, A3) right column: strong helicity and B) Right panel: Bq = WOG with Bl, B2, B3 the 
same as in A. Larger order modes were discarded. 






The magnetic contribution is: 
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TABLE II: First Panel: Intensity relationship between the tangential and normal to the field 
components of the eighteen first periods for Bq = WG and for weak (first column), moderate 
(second column) and strong helicity (third column) cases. The (P) indicates in phase and (IP) 
indicates inverted phase. Second Panel: Intensity relationship between the cylindrical components 
of the eighteen first periods for Bq = WG and for weak (first column), moderate (second column) 
and strong helicity (third column) cases. 
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TABLE III: First Panel - Stable case: Variation of the Radius with the Twist for weak hehcity 
b = 0.05 and Bq = lOG. Second Panel - Unstable case: Variation of the Radius with the Twist for 
moderate hehcity b = 0.5 and Bq = lOG. 

The heat conduction term results: 
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The new compressional contribution is expressed as: 
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were used. All the quantities were defined as in Paper I. 
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